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Abstract 

The structure of the Cauchy Horizon singularity of a black hole formed in a 
generic collapse is studied by means of a renormalization group equation for 
quantum gravity. It is shown that during the early evolution of the Cauchy 
Horizon the increase of the mass function is damped when quantum fluctua- 
tions of the metric are taken into account. 
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I. INTRODUCTION 



Recently much progress has been made in understanding the formation of singularities 
in realistic black holes. After the seminal work by Poisson and Israel the outcome of 
several investigations with spherical models (see for a general overview) was that the 
spacetime develops a null scalar singularity at the Cauchy Horizon (CH) whose subsequent 
evolution eventually stops at the final spacehke singularity at r = 

In particular the Petrov type D component \l/2 of the Weyl curvature diverges expo- 
nentially with advanced time at this lightlike hyper surf ace, although the "measured" tidal 
distortion is bounded. The metric tensor is regular in a suitable local chart adapted to the 
inner horizon and the metric perturbations are small. This scenario is likely to be essen- 
tially the same in more general contexts than spherical symmetry but the structure of 
the singularity is more complicated then since the square of the Weyl tensor C^^rxC^^'^^ is 
dominated by the radiative component \l/o^4 of a Petrov type N curvature 0. 

It is an interesting question whether quantum effects can modify the classical evolution 
of the fields in a significant way. As far as the classical evolution is concerned, causality does 
not permit our ignorance about the correct form of the djTiamics in the inner, Planckian 
curvature regions of the interior to infect the description of the overlaying layers in terms 
of classical general relativity. The radial coordinate r is in fact timelike in the interior of a 
spherical hole. 

This picture changes in quantum field theory because in loop calculations even states 
localized outside the light-cone have an impact on the value of the renormalized quantities. 
One can then imagine that the metric fluctuations near the inner horizon modify the infrared 
region where the Weyl curvature is still growing but it has not yet reached Planckian levels. 
In particular it is interesting to see if the presence of some "self-regulator" mechanism 
could prevent the local curvature from diverging at the CH. An indication has been given 
in 1^ where it has been noticed that the classical divergence of the mass function in an 
evaporating black hole (BH) can be damped out by the contribution of the blueshifted 
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influx of tlie Hawking radiation at late advanced times. More complete investigations in 
four dimensions ^ have been performed in the semiclassical approximation by considering 
a massless minimally coupled scalar field, but they were inconclusive about the "sign" of 
the quantum correction, i.e. about whether it would lead to a stronger or to a weaker 
divergence. 

The objective of the present investigation is to show that quantum fluctuations of the 
gravitational field indeed weaken the strength of the singularity at the inner horizon. This 
result is obtained by studying the running of the Newton constant at large momenta by 
means of the non-perturbative renormalization-group equation [^-[T0[| which governs the scale 
dependence of the effective average action for gravity ||lT[. is a Wilson- type effective 
action with a built-in infrared (IR) cut-off at the mass scale k. The functional is obtained 
by integrating out the quantum fluctuations with momenta between a fixed ultraviolet (UV) 
cutoff kuY and the variable IR cutoff k. In this framework, a renormalizable theory with the 
classical action 5* is quantized by solving the flow equation subject to the initial condition 
Tfcuv — ^ letting then fcuv oo, k ^ (after suitable renormalization) . 

What makes the effective average action an ideal tool for studying quantum gravity 
is the fact that this method can also be used in order to renormalization group-evolve 
(coarse grain) the actions of non-renormalizable effective field theories. In this case one 
assumes that there is some fundamental theory which has been "partially quantised", i.e. 
its quantum fluctuations with momenta from infinity down to a fixed scale k^y have been 
integrated out already. This leads to an effective action S^s which, when evaluated in tree 
approximation, correctly describes all phenomena with typical momenta of the order fcuv- If 
we are interested in processes at smaller momenta k < k\jy we can construct a new effective 
action, appropriate for the lower scale, by setting T^^^ = Scs and solving the flow equation 
for Ffc with this initial condition. It is clear that for effective theories the limit fcuv oo 
should not be performed; hence the non-renormalizability of a theory does not pose any 
problems in this context. 

Quite generally, the effective action F or the average action F^ encapsulates all physical 
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effects of a given theory. Once we have identified its leading terms for a given range of 
momenta, no other quantum corrections beyond those which are already contained in the 
running coupling constants parametrising the approximate form of Tk need to be taken into 
account. In the case at hand, a truncated derivative expansion in powers of the curvature 
tensor and its covariant derivatives is a sensible approximation as long as k « nip since 
the terms omitted are suppressed by higher powers of k/rup. (mp denotes the Planck mass.) 
Within this approximation, the most relevant term in the action is the Einstein-Hilbert term 
since it has the smallest canonical dimension. As a consequence, the most important effects 
of quantum gravity are encoded in the running of the associated coupling, i.e. Newton's 
constant. Thus, when we increase k from the classical (i.e, IR) domain to larger values in 
order to explore gravity at smaller distances, the first sign of a non-classical behaviour is 
a changing value of Newton's constant. In the present paper we investigate the regime of 
momenta where, on the one hand, the first quantum gravitational effects appear already 
while on the other hand higher order invariants (i?^ terms, etc.) are not yet important. 

In the following we shall consider Einstein gravity as an effective field theory and we 
identify the standard Einstein-Hilbert action with the average action r^^^^. Here /Cobs is 
some typical "observational scale" at which the classical tests of general relativity have 
confirmed the Einstein-Hilbert action. We assume that also for k > kohs, i-^- at higher 
energies, is well approximated by an action of the Einstein-Hilbert form as long as A; is 
not too close to the Planck scale. The two parameters in this action, Newton's constant and 
the cosmological constant, will depend on k, however, and the flow equation will tell us how 
the running Newton constant G{k) and the running cosmological constant A{k) depend on 
the cutoff. Their experimentally observed values are G{kobs) — Gobs and A(A;obs) = ^ohs — 0. 
The Newton constant defines the Planck mass according to rrip = 1/ V327rGobs- 

We fix a scale kuy ^ kohs in such a way that it is still sufficiently below the Planck 
scale so that Ffc is not yet very different from the Einstein-Hilbert action, but is already 
large enough for quantum gravitational effects to play a role. Wc start the renormalization 
group evolution at the scale kjjy with bare parameters G(A;uv) — G and A.{k\j\) — A which 
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should be thought of as functions of Gobs and Aobs- We shall use our result for the function 
G{k), k G [kohs, k\jy], in order to study the impact of the scale dependent Newton constant 
on the mass- inflation scenario. 

In a sense, we shall "renormalization group improve" the classical metric describing the 
late-time behavior of the spacetime near the CH. Our method is similar to the following 
renormalization-group based derivation of the Uehling correction to the Coulomb potential 
in massless QED ||T2[. One starts from the classical potential energy Ki(r) = e'^/Anr and 
replaces by the running gauge coupling in the one-loop approximation: 

e\k) =e\ko)[l-b\n{k/ko)]-\ h = e\k^)/Q'K\ (1) 

Hereby one may identify the renormalization point k with the inverse of the distance r 
because in the massless theory this is the only relevant scale. The result of this substitution 
reads 

V{r) = -e2(r^^)[l + 61n(ro/r) + 0(e^)]/47rr (2) 

where the IR reference scale tq = l/yUo has to be kept finite in the massless theory. Note 
that eq.(p is the correct (one-loop, massless) Uehling potential which is usually derived 
by standard perturbative methods |jl2|. Obviously the position dependent renormalization 
group improvement e^(A;), k oc 1/r encapsulates the most important effects which 

the quantum fluctuations have on the electric field produced by a point charge. We argue 
that analogous substitution G G{k) with an appropriate k = k{x^) yields the leading 
modification of the spacetime metric. 



II. QUANTUM GRAVITY EFFECTS BEHIND THE INNER POTENTIAL 

BARRIER 

Considering a spherically symmetric, charged BH the metric can be conveniently ex- 
pressed using the coordinates (a, 6 = 0, 1) on the radial two-spaces {6, 0) = const and the 
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function r(a;") that measures the area of those two-spheres whose hne element is r'^dD?. The 
metric element is then ds^ = gahdx"'dx^ + r'^dVl^. By defining the scalar fields f {x°') , m{x°') 
and —2k{x"') = df /dr through / = 1 — 2Gohs^/r + GohsG'^/r'^ the Einstein equations reduce 
to the two-dimensionally covariant equations 

r;ab + f^Qab = -4:nGohsriTab - QabT) 

R-2drK = SvrGobs (T - 2P) (3) 

where the static electro- magnetic field is generated by a charge of strength e and Tab is the 
stress-energy tensor of the matter field whose two-dimensional trace is T and tangential 
pressure is P. From the conservation laws one finds the following two-dimensional wave 
equation for the mass function 

□m = -leyrV^^GobsTabT"^ + 87rGobs/(P - T) 

+A7rr^G,^,KT - Airr^G^^^r^aT''^. (4) 

This latter equation is the key to understanding the phenomenon of the mass-inflation. The 
late time behavior of the external gravitational field produced during the collapse of a star 
is that of a (Kerr-Newman, in general) black hole of external mass mg perturbed by a tail 
of gravitational waves whose flux decays as ~ with p = 4{l + 1) for a multi-pole of 
order I. As a consequence of the boundary conditions set at the event horizon, the TabT"'^ 
interaction term between the influx and out-flux of gravitational waves scattered from the 
inner potential barrier triggers a divergent source term for the local mass function m{u,v). 
The outflow can be modelled as a radial stream of light-like material particles because of 
the infinite blue-shift near the Cauchy Horizon. It is possible to show |T|] that near the CH 

m{v,r) v'^e^"" {v ^ oo) (5) 

where v is the standard advanced time Eddington-Finkelstein coordinate, (kq denotes the 
surface gravity of the Reissner-Nordstrom static black hole that characterises the external 
field configuration.) 
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It must be observed that in Eq.(|^) the strength of the gravitational interactions between 
out-flux and influx is proportional to the Newton constant. Hence small changes in G due 
to renormalization effects are then exponentially amplified by the mass function like in a 
magnifying lens! In particular if gravity is asymptotically free the classical divergence of the 
mass function can be weakened by the decreasing of the Newton constant at small distances. 




FIG. 1. Penrose conformal diagram of a collapsing star. Note that the point H is not part of 
the manifold, but a singular point of this mapping. 

In order to discuss this phenomenon in the mass-inflation scenario we consider the model 
analysed in ||^ for the scalar field collapse although our result should not depend on this 
particular framework. We are interested in the asymptotic portion of the spacetime at late 
retarded times (the "corner" region near the point H in Fig.(|I|)) before the strong focusing 
region where r — > 0. The null Kruskal coordinates [/, V are thus introduced, being 

i^qU = — exp(/to'u), kqV = — exp(/tof) (6) 

were [u, v) are the retarded and advanced time coordinates. In a neighbourood of {U = 
— oo, V = 0) a.n approximate analytical solution of the Einstein equations and the wave 
equation for a massless minimally coupled scalar field $ can be found 0]. The explicit 
asymptotic expression for the metric is 
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ds'^ = -2—dUdV + r'^dQ'^ (7) 
r 
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r',-2Gobs[A{U)+B{V)] (8) 

where tq is the location of the CH in the static BH spacetime configuration. The dimension- 
less functions A{U) and B{V) are regular at the CH, A{—oo) = -B(O) = 0, but B diverges 
like l/V{—\n{—KoV)Y^^'^^ as ^ 0^ while A is positive definite and A is bounded. Even 
though the metric coefficients and the scalar field $ are both regular at the CH, the mass 
function is divergent for 0^ being 

m{U, V) ~ ^AB (9) 

We consider the evolution of the above geometry in the mass-infiating regime starting 
from a value of the coordinate V = Vin for which the mass function is already exponentially 
growing m{U, V)/mQ ^ 1 (we assume Tq ^ ~ rriQ ^ 1 in Planck units) but the curvature has 
not yet reached Planckian values. 

At this point we need an explicit expression for the running Newton constant. We use the 
result obtained in |TI[] where, for pure gravity, the evolution of has been obtained in the 



"Einstein-Hilbert approximation" where only the y/—g and y/—gR operators are considered 
in the renormalization group fiow. This amounts to truncating the space of all the actions 
to those of the form 

Tk[g,g] = (167rG(A;))-i| d^'x^i-Rig) + 2Aik)} + S,i[g,g] (10) 

where Sgf is the classical background gauge fixing term. For this truncation the fiow equation 
reads 

kd,T,[grg] = iTr[(rf [^?,^?] +7^r[^])-lA:9,7^r[^7] 
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-Tr 



M[g,-g]+nt[g])-'kd,nt[g]\ (ii) 



— (2) 

where (yf^i,, and A4 denote the background metric, the Hessian of F^ with respect to the 
"ordinary" metric argument g^^, and the Faddeev-Popov ghost operator, respectively. The 
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operators TZf^^ and TZf^ are the IR cutoffs in the graviton and the ghost sector, respectively. 
They are defined in terms of an arbitrary smooth function TZk{p^) (interpolating between 
zero for ^ oo and a constant oc k"^ at = 0) by replacing with the graviton and ghost 
kinetic operator, respectively. Inside loops, they suppress the contribution from modes with 
covariant momenta p < k. 

Upon projecting the renormalization group flow on the two dimensional space spanned 
by the operators y/—g and ^/—gR the functional flow equation becomes two ordinary dif- 
ferential equations for G{k) and A(A;). The equation for the scale-derivative of the running 
dimensionless Newton constant g{k) = k'^G{k) is found to be 

kdkgik) = [2 + r]ikMk) (12) 



where ri{k) = gBi/ (1 — gB2) is an anomalous dimension involving two known functions |TT[ 
of the cosmological constant, Bi and i?2, which depend on the choice for TZ^ip'^). Contrary to 
the running of the dimensionless gauge coupling e{k) in QED, the beta-function describing 
the running of the Newton constant is not universal. It is scheme dependent even in the 
lowest order of the loop expansion. In our framework this is reflected by the 7?.fc-dependence 
of rj. To lowest order of an expansion in powers of k/rrip one may ignore the impact of the 
running cosmological constant on ri{k) and set A(A;) ~ 0. Thus, returning to physical units 
and retaining only the leading term of the /c/mp-expansion the solution to eq.(|l^) reads 

G{k) = Gobs[l - wGobsfc' + O(fcVmJ)]. (13) 
For pure gravity one obtains u = ujg = —Bi{A = 0)/2 > which assumes the numerical 



value = 4(1 — tt/1AA)/tc for a standard exponential cutoff [|TT|. While u depends on the 
shape of TZ^, it can be shown that uj is positive for any choice of this function. Consequently 
pure gravity is "antiscreening" : Newton's constant decreases as k increases, i.e. it is large 
in the IR and becomes smaller in the UV. 

Eq.([TB|) is believed to be reliable as long as fcuv is still below the Planck mass. If 
k < k\jY << rrip the effect, in the renormalised system, of higher curvature invariants 
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such as i?^, R^pR^" or R^ which were omitted from the ansatz ([lOD is indeed smalL In 
fact, those invariants have been classified accordingly to their anomalous scaling dimension 
which characterise the linearised renormalization group flow near the Gaussian fixed point 
||1CI|| . The result is that the flow in the UV region is determined only by the "relevant" 
operators ^ and ^/gR, and that any other invariant with a higher canonical dimension is 



suppressed by additional powers of k/rrip |T^. We shall then assume that 



k < k\jY = nip /a (14) 

with a a fixed number well above unity. This defines the domain of validity of our approxi- 
mation. 

It is straightforward to include matter fields. In our model it might appear natural to 
keep the electro-magnetic field classical but quantise the full scalar field. The only effect on 



the running of G is to shift the parameter u. Using the same cut-off as above, one finds |T3 
uj = ujgs = 4/7r — 37r/72 which, again, is positive and leads to the same qualitative features 
as pure gravity. 

The running of G has dramatic consequences for the mass-inflation scenario. The leading 
quantum correction of the metric is obtained by replacing Gobs in eq. (|^) for the mass function 
by the running Newton constant G{k) with an appropriately chosen scale k. Since G{k) < 
Gobs for any value of A; > kohs "we conclude that the quantum corrections tend to damp the 
increase of the mass function. This qualitative conclusion is independent of the precise 
definition of the cut-off k. It is a rather robust results therefore. 



III. AN IMPROVED MODEL 

We are now going to implement this mechanism in an iterative calculation, where the 
zeroth-order solution for the metric is substituted into the running of G in order to calculate 
the first-order correction of the metric. For the sake of simplicity let us now consider the 
simpler case of the cross-flow model discussed in 
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The first question to be answered is what is the analogue of the identification k cx 
1/r which we used in QED. We are looking for a x'^-dependent cut-off k = k{x^) which 
respects general coordinate invariance and which measure the typical mass scale set by the 
curvature of spacetime. Since in the case at hand the metric is spherically symmetric, the 
natural candidate for the cut-off is the "coulombian" component of the Weyl curvature 
0: k^ oc |\E'2| = Gobs'm'{U,V)/rQ. More precisely we use the classical metric (zeroth-order 
approximation) to define the position dependent IR cutoff by 

k\V) = max{&^|^2|} = max{b'^Gohsm{U, V)/rl} (15) 

with b another fixed number much larger than unity and the maximum is performed over the 
region near U — > — oo. Here we are invoking a kind of adiabatic approximation where the 
use of a position dependent cutoff is justified because the mass function m{U, V) is almost 
constant on the length scales at which the eigenmodes integrated out are varying. (A similar 
approximation has already been used in [[^ in a semi-classical calculation.) From k{V) one 
obtains a running Newton constant as a function of the V coordinate 

G{V) = G,,,,[l-AB{V)] (16) 

where 

A = ub'max{Gl^J{U)/r',}. (17) 

It is now possible to evolve the classical geometry in eq.(|^ by considering the running 
Newton's constant in the Einstein equations. Within our approximation, the improvement 
amounts to replacing 

GohsTab = GohsT^{V) + GohsT°l^^{U) (18) 

with 

TjP = G(\/)T^^(K) + GobsTr (f/)- (19) 
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This modified energy-momentum tensor is then covariantly conserved since it satisfies 

(Trv)^^ = o. 

From the Bianchi identities one obtains the following wave equation for the mass function 
□ (Gm) = -167rV=^GXfeT'^' + e^f^V". (20) 

The general solution is uniquely determined once the value of the fields along the charac- 
teristic U = UiR and V = Vir are given. Asymptotically the improved metric is still of the 
form d^) but eq.(8) is now replaced by 

= rl - 2GiV)BiV) - 2G,^AiU). (21) 

By noticing that G approaches its bare value very rapidly one finds that the leading term (as 
V 0^) on the right hand side of eq.(pO|) is now given by the classically divergent TabT""^ 
contribution. Thus, after the inner potential barrier, one finds the solution 

m(f/, V)-^{1- AB)AB - ^^ABA (22) 

which replaces the classical expression (|^). Inserting the function B one has more explicitly 
for the leading term 

m{UV)--^(l - r23l 

^ ' ro ^ y(_ln(-ft:oV))(p+2)^y(-ln(-fi:o\/))(P+2) ^^''> 

This is our main result. It confirms our earlier conclusion about the damped increase of 
the mass function within an improved approximation which takes the back-reaction of the 
metric into account. 

A priori it might have appeared equally possible to perform the substitution Gobs ~^ G{k) 
directly in the Einstein equations. The identification k"^ = G{k)m/r^ leads to a non-linear 
equation for k that up to 0{k'^/rn^) is equivalent to k"^ = Gobs^^/^o which was used before. 
However it is important to observe that in this case a decreasing G leads to the additional 
effect of lowering the value of the surface gravity at the inner horizon which, too, damps the 
increase of the mass function. 
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It should also be stressed that the above results were obtained by integrating out only 
the field modes with momenta between h^J\^\ and m^/a. While lowering the IR cutoff even 
further is difficult from the technical point of view (the adiabatic approximation is not avail- 
able any longer) the monotonicity of the function in eq.(|T3|) suggests that taking additional 
modes into account will lead to an even stronger damping of the classical increase of the 
mass function. On the other hand, by adding further matter fields the antiscreening nature 
of the gravitational interaction could be destroyed in principle. (In ref. fl^ a condition on 
the number of the various species of fields implying > can be read off.) We nevertheless 
believe that the qualitative features of our discussion will hold for arbitrary matter systems 
with > 0. In particular for the matter system consisting of a massless minimally coupled 
scalar field considered in this investigation the quantum-corrected geometry is less singular 
than its classical counterpart. 



IV. CONCLUSIONS 



We have discussed a possible physical mechanism which has the effect of damping the 
classical increase of the mass function behind the potential barrier inside a realistic black 
hole. We believe that this mechanism is operative already for black holes with masses M » 
rrip and before the Cauchy horizon singularity is reached, i.e. in a regime of sub-planckian 
curvatures where it can be calculated reliably. On the basis of the present investigation 
we cannot make any claim about the fate of the singularity at the CH. However, one can 
speculate that if the decrease of G continues and that G ^ at the CH, the geometry of 
the spacetime near the late-time portion of the CH is regular with a sub-planckian Weyl 
curvature. In order to settle this issue completely a more complete calculation is needed 
and we hope to address this problem in the future. 

It is a pleasure to thank W. Israel for stimulating discussions. One of us (M.R.) would 
like to thank the Department of Physics of Catania University for the cordial hospitality 
extended to him while this work was in progress. He is also grateful to INFN, Sezione di 

13 



Catania, for the financial support which made this visit possible. 



14 



REFERENCES 
[1] E. Poisson and W. Israel, Phys. Rev. D 41, 1796, (1990). 

[2] Internal structure of black holes and spacetime singularities, volume XIII of the Annals 
of the Israel Physical Society, edited by L. M. Burko and A. Ori (Institute of Physics, 
Bristol, 1997). 

[3] A. Bonanno, S. Droz, W. Israel and S. M. Morsink, Proc.R.Soc.Lond. A 450, 553-567, 
London (1995). 

[4] A.Ori, Phys. Rev. Lett. 67, 1991, 789; A. Bonanno Phys. Rev. D 53, 7373 (1996). 

[5] P. R. Brady, S. Droz, S. M. Morsink, Phys. Rev. D 58: 080434 (1998). 

[6] R. Balbinot, P. R. Brady, W. Israel and E. Poisson, Physics Lett. A 161, 223 (1991). 

[7] W. G. Anderson, P. R. Brady, W. Israel, S. M. Morsink, Phys. Rev. Lett. 70, 1041 
(1993); R. Balbinot, E. Poisson, Phys. Rev. Lett. 70, 13 (1993); 

For an introduction see: M. Renter in Proceedings of the 5th Hellenic School and Work- 
shop on Elementary Particle Physics, Corfu, Greece, 1995 and |hep-th/9602012. 



[9] M. Renter, C. Wetterich, Nucl. Phys. B 417, 181, 1994; Nucl. Phys. B 427, 291,1994 . 

[10] A.Bonanno Phys. Rev. D 52, 969, 1995; A. Bonanno D. Zappala Phys. Rev. D 55, 6135 
(1997) and |hep-ph/961 1271 . 



[11] M. Renter, Phys. Rev. D 57, 971, 1998, and |hep-th/96050"30 



[12] W.Dittrich, M. Renter, Effective Lagrangians in Quantum Electrodynamics, Springer, 
Berlin, 1985. 

[13] L. N. Granda, Europhys. Lett. 42, 487, (1998) 
[14] D. Don, R. Percacci, |hep-th/97072"^ 



15 



